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Abstract 

We give a detailed presentation of a recently proposed mechanism of generat- 
ing the energy scale of inflation by loop effects in quantum cosmology. We discuss 
the quantum origin of the early inflationary Universe from the no-boundary and 
tunneling quantum states and present a universal effective action algorithm for 
the distribution function of chaotic inflationary cosmologies in both of these 
states. The energy scale of inflation is calculated by finding a sharp probability 
peak in this distribution function for a tunneling model driven by the inflaton field 
with large negative constant ^ of non-minimal interaction. The sub-Planckian 
parameters of this peak (the mean value of the corresponding Hubble constant 
H ~ 10~^mp, its quantum width IS.H/H ~ 10~^ and the number of inflationary 
e-foldings N > 60) are found to be in good correspondence with the observa- 
tional status of inflation theory, provided the coupling constants of the theory 
are constrained by a condition which is likely to be enforced by the (quasi) su- 
persymmetric nature of the sub-Planckian particle physics model. 

1. Introduction 

In this paper we give a detailed account of a recently proposed mechanism for 
generating the energy scale of the chaotic inflationary Universe by the loop part of the 
effective action in quantum cosmology 

Quantum cosmology became a theory of the quantum origin of inflationary Universe 
in early eighties due to the synthesis of the cosmological inflation with the idea of 
the quantum state 0, ^, generating the initial conditions for inflationary scenario. 
One of the main problems of this theory was the formulation of such a quantum state 
that could describe a very early quantum Universe, its evolution leading to the mod- 
ern observable large-scale structure. The inflation paradigm is very attractive because 



1 



at least heuristically it allows one to avoid applications of quantum gravity since the 
inflationary epoch has to take place at the energy scale or a characteristic value of the 
Hubble constant if = d/a ~ 10~^mp much below the Planck one mp = G^^"^. The 
predictions of the inflation theory essentially depend on this energy scale which must 
be chosen to provide a sufficient number of e-foldings N in the expansion law of a scale 
factor a{t) during the inflationary epoch, N = J^^ dt H > 60, and also generate the 
necessary level of density perturbations. This quantity, however, is a free parameter 
in the inflation theory, and, to the best of our knowledge, there are no convincing 
principles that could fix it without invoking the ideas of quantum gravity and cos- 
mology. These ideas imply that a quantum state of the Universe in the semiclassical 
regime gives rise to an ensemble of infiationary universes with different values of H, 
approximately evolving at later times according to classical equations of motion. This 
quantum state allows one to calculate the distribution function of this ensemble and 
interpret its probability maximum at certain value of H (if any) as generating the 
quantum scale of infiation. 

The implementation of this idea in the tree-level approximation of quantum cosmol- 
ogy P, 1^, |] has a controversial status and, in our opinion, is not satisfactory. The 
corresponding distribution functions are extremely flat |] for large values of H and 
unnormalizable at H ^ oo. This violates the validity of the semiclassical expansion 
underlying the inflation theory, since the contribution of over-Planckian energy scales 
is not suppressed to zero, and special assumptions are necessary to establish a Planck- 
ian "boundary" to protect semiclassical inflation physics from the nonperturbative 
realm of quantum gravity. Apart from this difficulty, the possible local maxima of 
the distribution function for the tree-level quantum states are either generating insuf- 
ficient amount of infiation violating the above bound [TU|, or generate too high level of 
quantum inhomogeneities and require unnaturally strong fine tuning (see Sect.2.1). 

The key to the solution of these problems, not resorting to the conjectures on a 
hypothetical over-Planckian phase of the theory, may consist in the semiclassical h- 
expansion and the search for mechanisms that could justify this expansion. Despite 
the perturbative nonrenormalizability of quantum gravity, this approach makes sense 
in problems with quantum states peaked at sub-Planckian energies. In particular, it 
would work in quantum cosmology with the no-boundary 0, ^ or tunnelling |^ wave- 
functions, provided they suppress the contribution of Planckian energies and generate 
the probability peaks at the lower (preferably GUT) scale compatible with microwave 
background observations. As shown in authors' papers |12], |l], |l3l, the loop effects 
can drastically change the predictions of the tree-level theory and really allow one to 
reach this goal. Moreover, as it was briefly announced in one can get a sharp 
probability peak in the distribution function of inflationary models with characteristic 
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parameters of GUT and, in this way, provide a numerically sound link between quan- 
tum cosmology, inflation theory and the particle physics of the early universe. Thus, 
the purpose of this paper is to give a detailed presentation of this work. 

The organization of the paper is as follows. Sec. 2 presents quantum cosmology as 
a theory of the quantum origin of the chaotic inflationary Universe. It gives a brief 
account of the quantum gravitational tunnelling underlying the no-boundary and tun- 
nelling wavefunctions of the Universe, discusses the model with nonminimally coupled 
inflaton fleld and presents a special algorithm for the one-loop distribution function 
of (quasi)-DeSitter models. Sects. 3-6 contain detailed calculations of various per- 
turbative contributions to this distribution function. We work within the double per- 
turbation theory: loop expansion in h up to the one-loop order and the expansion of 
the slow roll approximation up to the subleading order in the corresponding smallness 
parameter mp/|^|(^^ -C 1, where is a value of the inflaton scalar fleld and ^ = — 1^| is 
a big negative constant of its nonminimal coupling with curvature. Important differ- 
ence from our previous work where the distribution function was calculated only 
in the leading order of the slow roll expansion, is that now we flnd it in a subleading 
order. This leaves the main conclusions of [||] qualitatively the same and thus proves 
the stability of the leading order, although gives rise to certain quantitative correc- 
tions. Sects. 5 and 6 contain perturbative calculations of the one-loop effective action 
for generic set of flelds of various spins, contributing to the distribution function. In 
Sect. 7 we present a flnal answer for this function and flnd a corresponding probability 
peak that can be interpreted as generating the energy scale of inflation. In Sect. 8 
this result is used for the derivation of the selection criterion for viable particle physics 
models in the early Universe with nonminimal inflaton scalar fleld, apparently suggest- 
ing their (quasi)supersymmetric nature. This conclusion is based on the observation 
that the energy scale of inflation is suppressed relative to the Planck scale by the same 
small factor ~ 10~^ that determines a recently observed magnitude of the microwave 
background radiation anisotropy , provided that the two special combinations of 

coupling constants of the system satisfy certain restrictions. Sect. 9 contains concluding 
remarks. 



2. Quantum cosmology — the theory of quantum 
origin of the early inflationary Universe 

2.1. The no-boundary and tunneling quantum states 

It is widely recognized now that the inflationary scenario is one of the most promis- 
ing pictures of the early Universe 0. It can be described by the DeSitter or quasi- 
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DeSitter spacetime generated by an effective cosmological constant A, which in its turn 
is being generated by other slowly varying fields. Thus, in the model of chaotic inflation 
with the scalar inflaton fleld 0, minimally coupled to the metric tensor G^j^^ 

LiG,.,^) = G'/' I^^RiG,.) - i(V0)^ - f/(0)| , (2.1) 

the effective cosmological constant is generated by the potential of the inflaton fleld 
U{(f)). For this system rolling down from the potential barrier U{(j)) (which is supposed 
to be monotonically growing with 0) there exists the so-called slow roll approximation, 
when the non-stationarity of is much less than the rate of change of the cosmological 
scale factor a measured by the Hubble constant H = a/a. In this approximation 
equations of motion take the form 

^ ^^<if0, (2.2) 



H = H{<P) 



3H i 

\TTU{(f)) 



\ 3ml 



(2.3) 



so that an effective cosmological constant A = 3H^ is determined by the inflaton fleld 
potential. This potential is approximately constant during the inflationary stage due to 
the slow change of and only at the end of this stage decreases close to zero when the 
effective cosmological constant "decays" into inflaton oscillations, their energy being 
spent for the reheating of the Universe and its transition to radiation-dominated and 
then matter-dominated stages. The duration of the inflation stage usually measured 
by the number of e-foldings between the beginning tj and the end tp of inflation 

N= f^^dtH (2.4) 



ti 

determines the coeflicient of inflationary expansion of the model exp N and depends 
on the initial value of the inflaton fleld 0/. This dependence can be approximately 
obtained by changing the integration variable here to and integrating from 0/ to 
zero. This leads to a fundamental bound on the choice of inflationary model and initial 
value of inflaton N > QO M 



AT 0,) ^ ^ 
m. 



An ft^i H{(l)) 



dH{(P)/d(p 



> 60. (2.5) 



The role of quantum cosmology consists in the formulation of the quantum initial 
data for such a picture in the form of a particular quantum state of the Universe - 
the wavefunction ^{q) usually deflned on superspace of 3- metric coefficients and all 
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matter fields. The implementation of this idea was proposed in the pioneering works 
of Hartle, Hawking and Vilenkin ||, who suggested that such initial data (and 
the wavefunction ^{q)) correspond to the quantum gravitational tunneling that can 
semiclassically be described by the transition with changing spacetime signature. 

In the context of spatially closed cosmology the Lorentzian DeSitter spacetime 
can be regarded as a result of quantum tunneling from the classically forbidden state 
described by Euclidean DeSitter geometry. A simple picture of the tunneling geometry 
demonstrating such a mechanism is presented on Fig.l. The DeSitter solution of the 
Einstein equations with the cosmological constant A = 3H^ 

dsl = -dt^ + al{t) Cab dx^dx'', (2.6) 
aL{t) = ^ cosh (Ht) (2.7) 

describes the expansion of a spacelike spherical hypersurface with a metric of a 3- 
sphere a|(t) Cab with the radius (scale factor) aLit). Its Euclidean counterpart with the 
DeSitter positive signature metric 

ds"^ = dr^ + a^{T) Cab dx^dx^, (2.8) 
a (r) = 1 sin (ifr), (2.9) 

represents a geometry of a 4-dimensional sphere of radius R = 1/H with 3-dimensional 
sections (3-spheres) parameterized by a latitude angle 6 = Ht. These two metrics 
are related by the analytic continuation into a complex plane of the Euclidean time r 

mm 



T = n/2H + it, aL{t) = a{n/2H + it). (2.10) 

This analytic continuation can be interpreted as a quantum nucleation of the Lorentzian 
spacetime from the Euclidean one and shown on Fig.l as a matching of the two mani- 



folds ( p.6| ) - ( |2.9| ) across the equatorial section of the 4-sphere r = ti/2H (t = 0) - the 
bounce surface S^. 

Two known quantum states which semiclassically implement this mechanism, are 
represented by the no-boundary wave function of Hartle and Hawking ^ and the 
tunneling wave function known for historical reasons as a wavefunction of Vilenkin 
who pioneered the idea of quantum gravitational tunnelling in In the approximation 
of the two-dimensional minisuperspace consisting of the scale factor a and infiaton 
scalar field 

= («, 0), (2.11) 
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these wavef unctions '^nb{cl,<P) ^T(a,0) satisfy the minisuperspace version of the 
Wheeler-DeWitt equation and semiclassically represent its two hnear independent so- 
lutions 

*wij(a,0) ~e~^(^'<^), *T(a,0) ~e+^(^'<^), (2.12) 

where the Euchdean Hamilton- Jacobi function I{a, 0) of the model is calculated at 
a particular family of solutions of classical Euclidean equations of motion subject to 
special boundary conditions of Hartle and Hawking at a = (the "initial" point of 
the extremal) and boundary conditions (a, 0) at the end point of the extremal - an 
argument of the wavefunction. At a = the derivative of the scalar field with respect to 
the Euclidean time r should be zero while da/dr = 1 (r measures the proper distance), 
which is equivalent to the requirement of regularity of a 4-metric in the neighbourhood 
of the pole of a 4-sphere ( |2.8| ) - ( |2.9| ) at r = 0. In the leading order of the slow-roll 
approximation, when the inflaton field is constant, such a solution coincides with the 
exact round 4-metric ( p.8| ) - ( p.9| ) with the Hubble constant ( |2.3| ), and its Hamilton- 
Jacobi function equals 



^7rf/(0) 



lia, 0) = [l - (1 - if^(0) ay/'l , H\<^) = (2.13) 

When the point (a, 0) belongs to the region of minisuperspace below the curve (see 
Fig. 2) 

1 



the universe exists in the classically forbidden (underbarrier) state described by this 
Euclidean spacetime. Euclidean extremals, beginning at a = 0, have a caustic ( p.l4| ) 
and cannot penetrate into the region a > l/if(0) with real Euclidean time. However, 
they can be continued into this region by the analytic continuation into the complex 
time (|2.1CI| ) which generates the Lorentzian (imaginary) part of the Euclidean function 
J(a,0) 

J(a, 0) = J(0) ± iS{a, 0), a > l/i/(0), (2.15) 



^(«,0) = -^(^l0)«^-lj • (2.16) 

Here I{4>) is a Euclidean action of the theory with the Lagrangian ( |2.1D calculated on 
the gravitational half-instanton - the hemisphere (|2.8|) - (|2.9|) (0 < r < n/2H) 



^In the lowest order of the slow- roll approximation with constant (f) the problem is actually one- 
dimensional and Eq. ( 2.14 ) represents a set of turning points, however beyond this approximation 
this curve should, in fact, be replaced by the envelope of the family of Euclidean trajectories 
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This action determines the amphtude of wavefunctions ( p.l2| ) in a classically-allowed 
(Lorentzian) domain 

*7Vij(a, 0) ~ e~-^('^) cos (5(a, 0) + 7r/4) , (2.18) 
*r(a,0) ~ e+^('^) +''^(^''^), a>l/iJ(0), (2.19) 

which is interpreted in the tree-level approximation as a distribution function for 
the one-parameter ensemble of Lorentzian inflationary universes characterized by the 
Hamilton- Jacobi function ( p. 161) . The parameter enumerating the members of this en- 



semble is a value of the inflaton fleld or the corresponding Hubble constant H = H[(j)) 
and scalar curvature of the DeSitter space. Its quantum distributions for the no- 
boundary pNsi.'P) ^-iid tunnelling Pt(0) P quantum states read 

P^b(0) ~ e-2^(<^), pt{<P) ~ e+2^(<^). (2.20) 

The difference between these two wave functions and their quantum distributions 
consists in the different boundary conditions in superspace: while the tunneling state 
^T(a,0) at a > 1/H{(f)) contains only the outgoing wave and describes an expanding 
universe, the no-boundary wave function ^nb{0', 4>) in the Lorentzian regime represents 
the superposition of states of expanding and contracting cosmologies corresponding 
to the components of ( p.l8| ) of positive and negative frequencies with respect to the 
minisuperspace coordinate a. The tunneling wave function is deflned by the above 
mentioned outgoing wave conditions in the Lorentzian region of superspace and an 
additional condition of the ^-independence of ^T{a,(p) at a — > p], |1^]. For the 
no-boundary wavefunction there exists a more fundamental and model-independent 
prescription in the form of the path integral over regular Euclidean 4-geometries 

which in the tree-level approximation is dominated by the expression ( p. 181 ) - a 



contribution of the saddle point of this integral - the Euclidean-Lorentzian extremal 



m - m- 

The distribution functions Pnb{(P) and pr(0) describe the opposite results of the 
most probable underbarrier tunneling: to the minimum and maximum of the inflaton 
potential U{(f)) > correspondingly (though in the latter case a minimum U{(j)) = 
does not belong, strictly speaking, to the domain of applicability of the slow-roll 
approximation) . 

Equations given above apply to the model ( |2.1| ), however they can also be used in 
the theory with non-minimally coupled scalar inflaton 

L{g,.,v) = g'" [^R(9,.) -\iv'R{9,.) -\{^vf -\n\' -\^'\^ (2.21) 
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provided L{G^i,, 0) above is viewed as the Einstein frame of the Lagrangian L{ g^,,, ip ) 
with the fields = {{l + '&Ti\^\(f'^ /rn?p)g^y,(j){ip)) related to [g^y.ip) by the known 

conformal transformation of the metric and the reparametrization of the scalar field 



19], |21|. For a negative nonminimal coupling constant ^ = — 1^| this model easily 
generates the chaotic inflation scenario with the following inflaton potential in the 
Einstein frame parameterization 

•t'^'t'if) (l + 8i|{ Ij^Vm?,) 

including the case of symmetry breaking at scale v with m? = — Af ^ < in the Higgs 
potential A(v9^ — v^^/A . At large if the potential ( p.22|) approaches a constant and 
depending on the parameter 

.2 O^ICU,2 



S ^ = (2.23) 

Amp rrip 



has two types of behaviour at the intermediate values of the inflaton field. For 6 > —1 
it does not have local maxima and generates the slow-roll decrease of the scalar field 
from its initial value ipi leading to a standard scenario with a finite inflationary stage 
and approximate e-folding number 

^(„).(^)^^M±M. ,2.24) 

Knip/ 1 + 



For 5 < — 1 it has a local maximum at (p = m/ y X\l + 6\, including the case of zero 

A when ip = nip/ -y/87r|^|, and due to a negative slope of the potential leads to the 
inflation with infinite duration for all models with the scalar field growing from its 
initial value ipi > ip io infinity. 

The tree-level distribution functions ( p.20| ) for such a potential do not suppress the 
over-Planckian scales and are unnormalizable at large dip pi\ib,t{v) = oo, thus 

invalidating a semiclassical expansion. Only for A = the normalizability takes place in 
the tunnelling case with |^| 7^ 0, Pt(v^) ~ exp[— 487r^|^pv9^/m^], ip ^ 00, but this fine 
tuning is too strong and can hardly survive renormalization of A by quantum effects. 

Only for a tunneling case with 6 < —1 the distribution Pt{4') has a local peak at the 
maximum of the potential (|2.22| ) (p, which could serve as a source of the energy scale of 



inflation at reasonable sub-Planckian value of the Hubble constant. However, this peak 
requires a large positive mass of the inflaton field > Amp/(87r|,^|), which is too large 
for reasonable values of ^ = —2 x 10^, A = 0.05 |19[. Another (and maybe more serious) 
difficulty with this inflation scenario starting from the maximum of inflation potential 
is that according to (p.2|) it begins with zero ~ = and generates infinitely large 
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quantum inhomogeneities (inverse proportional to 0]) which are incompatible with 
the observable large scale structure of the Universe. All this makes questionable the 
attempts to arrange the quantum origin of our Universe at the tree-level theory and 
serves as a motivation for considering the loop effects. 

2.2. One-loop distribution function of the inflationary cos- 
mologies 

Note that the calculation of the tree-level distribution does not require the knowl- 
edge of the correct probabilistic inner product of cosmological wave functions. It is 
enough to calculate and square an amplitude of the wave function, which due to the 
peculiarities of the model is a function on the section of the two-dimensional minisu- 
perspace, transversal to the coordinate a, usually playing the role of time. Therefore 
the obtained distribution function is defined on the physical subspace of a correct di- 
mensionality - one- dimensional space of spatially homogeneous inflaton field. Beyond 
the tree-level approximation the situation changes: one needs the knowledge of the 
wave function with the preexponential factor in the needed approximation, knowledge 
of the correct inner product and the extension beyond the minisuperspace approxima- 
tion, because the distribution function contains now a non-trivial contribution from 
integration over inhomogeneous quantum fields frozen in the tree-level approximation. 
At the one-loop level which we shall study here it is enough to consider these fields in 
the linear approximation. Setting of the problem in the model of chaotic inflationary 
Universe consists in the minisuperspace model with the scale factor a and spatially ho- 
mogeneous scalar inflaton ip and with inhomogeneous fields of all possible spins treated 
as perturbations on this background. Together they form a superspace of variables 

q = (a, V2, V5(x), V^(x), A„(x), V'a(x), hab (x), ...). (2.25) 

On this superspace we shall have to calculate the no-boundary and tunnelling wave- 
functions ^nb,t{(i) and then, by using a proper physical inner product, calculate the 
distribution function of the collective variable ip. To make the latter step and even 
to give a rigorous definition of this distribution, it is better first to make a quantum 
reduction to the wave function of physical ADM variables ^, which simultaneously 
disentangles time t (initially parametrized among the superspace variables (|2.25|) p3[| ) 

q^i^,t), *(g)^v]>(e,t). (2.26) 

Strictly speaking this reduction is not consistent (globally on phase space of the 
theory), and a complete understanding and the interpretation of the cosmological wave- 
function might be reached only in the conceptually more advanced framework (third 
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quantization of gravity theory, refined algebraic quantization |2J, |2^, etc.). Althougli 
tliis framework still does not have a status of a well-established physical theory, there 
exists a good correspondence principle of this framework with the quantization in 
reduced phase space for systems with a wide class of special (positive-frequency) semi- 
classical quantum states. For these states the conserved current of the Wheeler-DeWitt 
equations perturbatively coincides with the inner product of the ADM quantization and 
thus can be used for the construction of the probability distribution (for a perturbative 
equivalence of the ADM and Dirac- Wheeler-DeWitt quantization of gravity for such 
physical states see p6|, |12|)- The tunneling wavefunction belongs to such a class 
of states, while the no-boundary one does not and should be supplied with additional 
(third quantization) principles to be interpreted in terms of the probability distribution 
of the above type. 

It is plausible to make this reduction separately in the minisuperspace sector of the 
full superspace (a, if) and its sector of inhomogeneous modes. We choose an inflaton 
field as a physical variable whose distribution function we will calculate, while the 
solution of classical equations of motion (|2.7| ) with H = H{(f) will be considered as a 
gauge 

X^(a, ip,t) = a- - J-^cosh(ff(¥;)t) = 0. (2.27) 

It simultaneously plays the role of the parameterization of minisuperspace coordinates 
in terms of the physical variable 0. The ADM reduction for linearized inhomogeneous 
modes of fields boils down to the choice of their transverse (T) and transverse-traceless 
components (TT), so that the full set of physical variables reads 

= /), / = (^(x), ^(x), Ar(x), ^r(x), hl^i^), ...). (2.28) 

At the quantum level the ADM reduction can be easily carried out by the method 
described in [|^, ^ for the tunneling state (|2.19|) . However it stumbles upon the 



problem of positive and negative frequency components for the Hartle-Hawking wave- 
function ( |2.18| ) and in the gauge ( |2.27| ) encounters the analogue of the Gribov copies 



problem, corresponding to these components. As discussed in p9[, these copies are 
an artifact of using inappropriate gauge, whose surface intersects twice the classical 
extremal ( |2.7| ) of one and the same Universe before and after its bounce against the 
minimal value of a cosmological radius a = 1/H{(p). This implies a dubious interpreta- 



tion of (|2.18| ) as a superposition of two simultaneously existing states of expanding and 



contracting Universe. This problem can be resolved at the fundamental level by the 



^This gauge is very convenient because it approximately corresponds to the choice of the proper 
time with the lapse function N-^ = 1 . 
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transition to quantization in conformal superspace in the framework of the York gauge 



29| , but this framework is not yet developed to be a workable technique. However, in 
the present model in the semiclassical approximation it is enough merely to consider 
the quantum ADM reduction for a separate positive frequency (or negative frequency) 
component of ( 2.18|) . 

Thus, semiclassically for both cosmological states the quantum ADM reduction 
boils down to obtaining the corresponding wave function of physical variables \l/(^, t) = 
"^{(f, f\t). Then, the distribution function of if should be regarded as a diagonal element 
of the density matrix of this pure state Trj|\E')(^|. It can be obtained from |\E') = 
f \ t) by averaging over the rest of the modes of physical fields / 



p(0|t)= df^*{<PJ\t)<ili^,f\t), 



(2.29) 

and does not reduce to a simple squaring of the wave function. 

The calculation of the one-loop no-boundary and tunnelling wavefunctions pertur- 
batively in inhomogeneous modes / on the Friedmann-Robertson- Walker background 
was carried out in many papers ^ |3^, 0, It can be based on the 

path integration over the fields regular on the Euclidean spacetime with metric (|2.8| ) 
- (|2.9| ) or by using the known one-loop approximation for the general solution of the 
Wheeler-DeWitt equations ||2^, Then both wavefunctions turn out to be Gaus- 

sian in the variables / - their Euclidean DeSitter invariant vacuum [|1^]. Therefore the 
integration over / in (^.29|) is trivial and leads to the fundamental algorithm which is 



valid for both no-boundary |TT|, [T^, |TB|, |I[ and tunneling |Q quantum states 

PNsA^lt) = e^^^^^^ " (2.30) 

It involves the doubled Euclidean action on the hemisphere with the metric ( p.8|) - 
(|2.9|), the linearized mode of the homogeneous inflaton field v^(t) - the basis function 
of the wave equation on the Lorentzian DeSitter background (|2] 



be obtained from the regular Euclidean linearized mode 







by the analytic continuation ( p. 101) 



Si 6^ 



0, 



{n/2H + it)]* 



which can 
(2.31) 

(2.32) 
(2.33) 



Here /[^] and S[i] are the Euclidean and Lorentzian action functional reduced to 
the ADM physical variables, so that S'^I[i]/SiSi and S'^S[i]/6C, 5^ are the operators of 
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their physical inverse propagators. The mode v^^it) is normahzed in ( p.30| ) to unity 
with respect to the Wronskian inner product in the space of solutions of this wave 
equation 

A^ = iv;{t)W v^{t), (2.34) 

* — > < — - 

where W=W — PI/ is a corresponding Wronskian operator linear in time derivative. 
The one-loop contribution is the same for both quantum states and is determined by 
the Euclidean effective action of all physical fields C,{x) 



ri_ioopM= ^Trln-'^'^t^- 



2 S^S^ 



(2.35) 

DS 



This effective action is calculated on the DeSitter instanton - 4-dimensional sphere of 
the radius 1/H{<f) - and, therefore, is a function of <y9 - the argument of the distribution 
function. Such a closed Euclidean manifold is obtained by the doubling of the half- 
instanton - two hemispheres match each other along the equatorial hypersurface 



(on which a quantum transition with the change of signature takes place). Graphical 
illustration of the procedure of calculating the distribution function is given on Fig. 3. 
The wave function and its conjugate, participating in the scalar product ( |2.29|) , can 



be represented by two Euclidean-Lorentzian manifolds. When calculating the inner 
product, due to implicit unitarity of the theory, the contributions of Lorentzian regions 
cancel each other and the result boils down to the Euclidean effective action calculated 
on a closed instanton obtained by gluing the two hemispheres of the above type [|T^, 

The above simple scheme holds only in the lowest order of the slow-roll approxima- 
tion, when the inflaton scalar field is constant on the solution of classical (Euclidean and 
Lorentzian) equations of motion and generates the effective Hubble and cosmological 
constants invariable in time. This is the case of the so-called real tunnelling geometry 



35| . In the chaotic inflation model, however, the inflaton field is varying with time 
which makes the spacetime geometry deviating from the exact DeSitter one, and corre- 
spondingly the analytic continuation from the Euclidean regime to the Lorentzian one 
makes the extremal complex. The general case of complex gravitational tunneling was 
considered in much detail in [|r^ and for this particular model in [^. The modifica- 
tion due to the imaginary part of the complex extremal looks as follows. The complex 
minisuperspace extremal Q{z) = {a{z),ip{z)) in the complex plane of the Euclidean 
time z = T + it should satisfy the following boundary- value problem 



6I[Q] ^ difiz) 



6Q{z) ' dz 



0, a{z) z + 0(^2), z^O, (2.36) 



=0 



Q{z+) =q+ = {a+, <^+) (2.37) 



12 



with the Hartle-Hawking no-boundary conditions at z = and real boundary data 
(a+, ip^) at the boundary of the spacetime ball z = z+. The Euclidean action = 
I[Q{z)] calculated on the solution of this problem is in general complex but cannot be 
decomposed as before in the purely real contribution of the Euclidean part of the full 
manifold and imaginary contribution of its real Lorentzian section. Thus, the algorithm 
(|2.30D still holds, but with I{(p) replaced by the real part of the complex action 



3?e I[Qiz)] = g?e I{a+, (2.38) 
and with all the other quantities calculated on the background of this extremal . 



The final real point of the extremal (a+, v?+) should be a subject of the ADM reduction 
which identifies (p^ with the physical variable and expresses a_|_ as a function of the 
physical time and <y9_|_ in some gauge, like ( p.27| ). In what follows we shall find 



such a complex extremal in the first subleading order of the slow-roll approximation 
and calculate the corresponding distribution function. The model we consider will 
be a chaotic inflationary cosmology with the infiaton field non-minimally coupled to 
curvature with a large negative coupling constant ^ = — As we shall see, a small 
parameter of the slow roll expansion in this model turns out to be inverse proportional 
to 1^1, "mp/l^lv^^ <^ 1. This choice is justified by the fact that this model with |^| ~ 
2 X 10'^ is regarded as a good candidate for the inflationary scenario compatible with 
the observational status of inflation theory ||19||. 



3. Non-minimal infiaton scalar field: tree-level ap- 
proximation 

3.1. Perturbation theory for Euclidean classical solutions 



Let us find the classical solution of the problem ( p.36|) in the model (|2.21|) with 



the nonminimally coupled infiaton field. We shall begin solving this problem on the 
Euclidean segment of complex time z = t. The Euclidean action of this model in the 
minisuperspace of homogeneous scale factor a, infiaton ip and lapse looks like 



-6{W^ - imV - jAv,"} . (3.1) 
where the dot denotes the differentiation with respect to proper Euclidean time d/Ndr. 
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The first order variational derivatives of this action with respect to and give 
the Euchdean equations of motion 



a 



1 



a 
a 



0, 



(3.2) 



I - 



N \ 2\ An 



— m 09 Xip 

2 ^ 4 ^ 



0, 



(3.3) 



where the normahzation factor A = 27r'^Na^. Eq. (|3.2|) is a dynamical equation of 
motion for the scalar field if, while Eq. (|3.3|) is a Lagrangian version of the Hamiltonian 
constraint. There is no need to write down the variation of the action with respect to 
the scale factor a because this equation is a consequence of ( p.2|) - (|3.3| ). 

We now develop the perturbation theory for the above equations starting with the 
DeSitter solution as a lowest order approximation of a constant scalar field ip. This 
solution in the cosmic time gauge = 1 has the form 



(0) 



ifo = const, a 



(0) 



sin[i^ (v9o) r] 



(3.4) 



where the bracketed superscript denotes the order of perturbation theory and ipo = ip{0) 
is the initial value of the infiaton scalar field at r = and a = 0. Substituting it into 
Eq. (^.31) we obtain H{(f) as a following function of the scalar field: 



mp{l 



25) 



+ 



m. 



(3.5) 



On the solution (|3.4| ) the equation (|37^ ) is satisfied exactly, while the equation (|3.2| ) 
holds only approximately 



T (0) 



(3.6) 
(3.7) 



with the effective mass in the dynamical equation for infiaton field M'^{ip) which is 
small in view of (|3.5[) (remember that = — < 0) in the limit of large |^| and 



(3.1 



provided the parameter 6 defined by Eq. (|2.23| ) is bounded (in the chaotic inflation 
model with nonminimal infiaton the constant |^| is usually chosen of the order of 
magnitude of ratio of the Planck scale to the GUT scale v = —m?/X, so that 5 <^ 1). 
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This property of cancellation of the leading in contributions to (|3.71 ) underlies the 
slow-roll approximation in this model of the nonminimal and nonlinear infiaton field. 
This approximation works for large values of the infiaton, its inverse playing the role 
of smallness parameter. As we shall now see, big negative ^ further improves this 
expansion which actually takes place in powers of w?p/\^\ip^. 

To find the first subleading order of this expansion, we have to expand the equations 
of motion (|3.2|) -(|3l3|) up to the first order in perturbations 



/,Ar<^(5v9 + /,Ara5a = -/,Ar = 0, (3.9) 

I,^^5^ + I,^a5a= (3.10) 

and solve this linear system for 5ip and 5a. Here we use an obvious notations for the 
second order variational derivatives of the action, which are the differential operators 
evaluated at the lowest-order solution ( |3.4| ). With the choice of a new (angular) variable 
on the DeSitter sphere, replacing the Euclidean time, 

e = H{^^)T (3.11) 

these operators take the form 

/w--^/^^{^ + 3cot.A_(^_24«)), (3.12) 

,,^„._^^^5ip{^^2cot4.3}.il,„ (3.13) 
= AH^ {65^„ (l - cot « A) + !E?!|±M} , (3.14) 

I.Na= -—{'iml/AT^ - Qi^l) il + cote—y (3.15) 

with the normalization factor A, Hubble constant H and scale factor a taken in the 
lowest-order approximation. 

For large ipo and in view of (|3.5| ) the last two equations can be simplified 

-^/,iv<^(5v? = -e^y^o^^^?^-^^^ + 0(mp/(/?o), (3.16) 
AH^ smu au cos 6 

— -J,^, Sa = -Q^^lH-^ — (1 + Oiml^D) 3.17 

AH"^ sm 6 aU cos ^ ' 



and used in the linearized constraint equation ( p.9| ) to give 
d 5a sin 9 d 5ip 



dO cos 9 ^oHq dO cos 9 



+ 0{m%/^l). (3.18) 



15 



This equation in its turn can be used in the dynamical hnearized equation (|3.1CI| ). 
The latter in view of (p.l2| )- (|3.13| ) can be rewritten in the form 

<^ — + 3 cot ^— + 24^ \ 6ip + Q^ifoH <^ — + 3 cot ^ - tan^ 



^de^ dO "J "'^^ {de J sin^d^cos^ 

= ^ + 0{mll^ll (3.19) 

which allows one to exclude 5a on account of (|3.18|) and thus arrive at the closed 
equation for 5({) 

Its solution satisfying the Hartle-Hawking boundary conditions 0, |[ {d/ d6)5ip{Q) = 
0, 6ip{0) = 0, reads (we use ( ^3)) and (|37^ ) 



7r(/?o 1 — 64\ 2 4 2y 

The corresponding perturbation of the scale factor 5a can then be obtained by inte- 
grating the equation (|3.18|) with zero initial condition at 6' = 

5a = — — — cos y tan m cos tan -tanc? 



■KiplH{^a) 1-6^ V 2 4 2 

e 1 6' 
'a^2 



tan-- Trf^'lncos- ) +0(m|/y?^). (3.22) 
2 JO 2 / 



The results of this section can be summarized in the one-parameter family of Eu- 
clidean extremals enumerated by the initial value of the infiaton field <^o 

Lp{T,ipQ) = ipQ + 5ip{0,ipQ), (3.23) 
«(^'<^o) = 777^sin^ + 5a(^,(^o), 9 = H{^o)t (3.24) 

with perturbations 5(f = 0{mp/ (po) and 5a = 0{m^/(p^) given above. These pertur- 
bations are small for large yjo- The corresponding smallness parameter is m?p/ip'^ (one 
should remember that for a bounded range of 9 and in the leading order ~ (/jq, a ~ 
V^o^). A direct inspection of these perturbations also shows that for large negative 
constant of nonminimal coupling |,^| ^ 1 the actual smallness parameter is mp/|^|(y9Q, 
and the final conclusions of this paper will be obtained for 'rn?p/\^\ipl <^ 1. 
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3.2. Superspace caustic and complex extremals 

Once we have a perturbative solution parametrized by the initial value of the 
scalar field ipo, we can now perturbatively solve the boundary value problem ( ^.37] ) 
which takes the form of two equations 

Lpo + 5ip{6+,ipo) = ip+, (3.25) 
— — -sme+ + 6a{e+,(po) = a+, (3.26) 



for {ipo, 6+) in terms of (</)+, ci+). Like (|3.23| )- (|3.24|) this solution can be obtained in a 
subleading approximation in the form 

y,o = ^r+^i^ e, = 9?^ + 9^^\ (3.27) 

where the superscript in brackets denotes the order of perturbation theory. 
The lowest order approximation 

<^[,°) = (3.28) 

= arcsin [if (</?+) a+] (3.29) 

shows that only the points (</)+, a+) lying below the curve 

^ (3.30) 



can be reached by Euclidean trajectories (with the no-boundary initial conditions) in 
real time. This curve consists of the points of maximal expansion of the Euclidean 
model at 6 = -njl and can be regarded as a caustic of the family of solutions in the 
lowest order of the slow roll expansion. The points of two-dimensional minisuperspace 
above this curve can be reached only in complex Euclidean time because 



vr 



Qf = -^^E{^^)U, (3.31) 
Htj^ = arccosh [H{(p^)a^] , (3.32) 

where t+ can be interpreted as the Lorentzian time in the DeSitter space nucleating at 
t_i_ = from the Euclidean DeSitter hemisphere parameterized by the angular coordi- 
nate 6 = H{ip^)T, < 6 < 7r/2. Thus in this case the combined Euclidean-Lorentzian 
evolution takes place on a contour in a complex plane of time, consisting of the Eu- 
clidean < r < n/2H{{p+) and Lorentzian r = 7r/2if (</?+) + it, < t < segments. 

In the lowest order of the slow roll approximation the point of Euclidean-Lorentzian 
transition coincides with the point of maximal expansion of a on the Euclidean tra- 
jectory where the both velocities if and a vanish. This situation corresponds to the 
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so-called real tunnelling |^ when the analytically matched Euclidean and Lorentzian 
extremals both have real values of the configuration space coordinates. Beyond the 
leading order of the slow roll approximation this property does not hold. 

The first sub leading approximation for {ipo, 9+) can be obtained by making in Eqs. 
(|3.25| )- (|3.26|) the first-order iteration in 6ip and 6a 

rf) = -M^f (3.33) 
'^'^ ^ (i/(¥.+)<5a(^f ,¥.+) +sin^f <5<^(^f ,y,+)/y,+) = 

1 0^°^ ro^"^ b\ 

tan^ + / ^ In cos- +0{m%/^X). (3.34) 
I I Jo 1 1 

In this approximation the point of maximal expansion (extremum of (|3.26| )) at 6m 
differs from tx/2 

TT 

Om = -+Sm, (3.35) 

"^'^ ' + ^ ^ ' ■ " "ln2 + GV (3.36) 



7r(^^ 1 - 6^ V 2 8 2 

where G is the Catalan constant {G = J^^^ dOlncotO = 0.915. . .). The set of these 
points, however, does not form the caustic curve: in contrast with the lowest order of 
the slow roll approximation, in which the dynamics is actually one-dimensional, 0o = 0, 
now the both minisuperspace coordinates depend on time and the equation that deter- 
mines the envelope of the 9?o-parameter family of extremals (|3.23|) - (|3.24|) is given by 



the following requirement. At the envelope curve the vector tangential to the extremal 
(0, d) and the vector transversal to their one parameter family {d(p/dipo,da/dipo) are 
coUinear, i.e. 

. da . d<^ 

(p- a-— = 0. (3.37) 

difo d(po 

From this equation it follows that the extremal starting at r = with the initial scalar 
field ipo hits the caustic at 

Oc{^o) = l + ec{vo), (3.38) 

mp 1 + 5/1 TT TT 



The caustic curve when parametrized by the value of ipo can be obtained by substituting 
9c{(po) in ( |3.23| ) and ( |3.24| ). Then, with ipo excluded in terms of (p from the second of 
the resulting equations, the first of them becomes the equation of the caustic as a 
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graph of a against (p. Simple calculations show that the first order corrections cancel 
out and the caustic equation remains the same as the leading-order one (|3.30| ) a = 
1/H{(p) + 0{mp/(p^). The meaning of this curve implied by the equation ( |3.37| ) is that 
its points belonging to a particular extremal a = a(r, yjo), f = ^ijiVo) are no longer 
in one-to-one correspondence with the value of the Euclidean time r and initial value 
of the scalar field ip^ enumerating the extremals 

777^ = 0- 3.40 

One cannot go beyond this curve in real Euclidean time: the continuation in real 
> 6'c(v?o) results in the trajectory bouncing back to the Euchdean domain under the 
caustic curve a < 1/ H{(f)+0{m'^/(p^). Similarly to the lowest order approximation the 
points beyond this curve can be reached only in complex time. Now, however, we have 
two important differences. Firstly, the both velocities (a, ip) at the nucleation point (as 
anywhere else on every extremal) are nonvanishing simultaneously which results in the 
complex valuedness of the extremal ending at the point (</?+, a+) beyond the caustic 
(the discussion of matching conditions at the Euclidean-Lorentzian transition for this 
case of complex tunnelling can be found in the previous authors' work |]T3|). And, 
secondly, the point of nucleation (the corresponding value of the Euclidean time) is a 
nontrivial function of the end point of the extremal. To show this take a+ > 1/H{(f^) 
and calculate the final value of the complex Euclidean time 6*+ given by Eqs. (|3.27|) , 
(|33lD , dpI), (|01) . It has a real part 

3fte e+ = ^jv(<^+, t+) = l + eN{v+, t+), (3.41) 

mp 1 + 6 f 1 IT 71 

= ^1^ [-2cosh[H{^^)U] + 8 " 2 2 + ^ 

1 /■*+ \ 

+-H{lp+) / dt arctan sinh [ H{ip+) t]] , (3.42) 

2 Jo J 

which can be identified with the nucleation point provided we choose the contour of 
complex "angular" time 6 joining the points 6 = and 6*+ with the union of two 
segments < 9 < 9^ and 9 = 9^ + iH{{p+)t, < t < t+. In contrast with the 
case of real tunnelling this is just a convention, because on both segments the fields 
are complex and analyticity does not give any preference to this particular choice. 
Comparison of Em, Sc and shows that all these deviations from the lowest order 
nucleation point 9 = 11/2 are different: Sm < and En = £c only at ^4. = 0, that is 
when the end point of the extremal lies on the caustic, while for positive Lorentzian 
time Sat > 0- 

''Note that our perturbation theory for 61^ and other Lorentzian quantities in its present form is 
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Once we have a solution for equations (|2.36|) - (|2.37|) we can calculate the corre 



sponding Hamilton- Jacobi function and, in particular, its real part (|2.38| ) Sfte /(a+, 93+). 
When calculated in the lowest order approximation on the solution ( |3.4| ) the doubled 
real part of the action obviously coincides with the Euclidean action on the full four- 
dimensional sphere of the radius 1/H{ip+) and constant scalar field ip = ip+. When 
expanded in powers of nip/ip'^ it equals 

2Ke /(a+, = Jo + A + o ('^] , (3.43) 



967r2|e|2 



A 



(3.44) 



A = JJ^amim, (3,45) 
A 

Since we want to have this quantity in the first sub leading approximation in mp/y)^, 
a priori we have to include corrections to this result linear in 6ip and 6a. By direct 
calculations these corrections can be shown to vanish due to certain intrinsic cancel- 
lations. The mechanism of these cancellations follows from the parametrized nature 
of the gravity theory and looks as follows. Using the notations of (|2.36| )- (|2.37| ) (with 
complex 6 replacing z) we can write down the total action on the extremal subject to 
boundary conditions g+ at 6*+ as an integral over 6 of the corresponding Lagrangian 

leAQi^)] = r deL{Q,Q) (3.46) 

•J 

with a subscript indicating the upper limit of integration ^+ over the complex angular 
"time" which is a solution of the boundary condition ( p.37|) . We have a solution of clas- 
sical equations of motion as a perturbation expansion (5(6*, v^o) = Q^^"* (^5 V^o) + V'o) 
and in two subsequent orders of this perturbation theory the boundary conditions have 
the form 

Q°(^f,^r) = g+ (3.47) 

Q(°)(^r + e^l\ + rf^) + mefJ^') = (3.48) 



whence it follows that 

= 0, (3.49) 



valid only for small values of the Lorentzian time t, H{ip^) i+ ~ 1, to be valid in the slow roll limit 
of big H{ip^). In the Euclidean context one of the motivations for introducing the coordinate 9 — Ht 
instead of r was the fact the range of 6 below the caustic is bounded by 7r/2 which is different from 
the Lorentzian domain where the hyperbolic "angle" Ht is unbounded from above. 
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where the dot denotes the derivative with respect to 9. Now, the total action ( p.46| ) 
calculated in a first subleading order of perturbation theory reads 



Q 



[0)1 



n(0) 



3(1) 



dL_ / ggw 

[dQ \ 



6Q 



(3.50) 



(here the first order variation of the functional argument of the action reduces to the 
standard surface term at because this variation is being calculated at the solution 
of the classical equations). Then from ( p.49|) the total first order correction induced by 
SQ reduces to 



L 








(3.51) 



and vanishes because the coefficient of here in a reparametrization invariant theory 
boils down to a Hamiltonian constraint identically satisfied for a classical background. 
Thus the final expression for the real part of the classical action reduces in the sub- 
leading approximation to ( p.43|) -( p^45|) . 

The absence of the first-order corrections in ^el{a+,ip+) due to 6Q guarantees 
the tree level unitarity of the theory - the time independence of the semiclassical 
wavefunction amplitude. As it follows from (|3.43| ) it depends only on the final value 
of the scalar field while the contribution of 6Q could have introduced a nontrivial 
dependence on a+ (or t+). This property, obvious for real tunnelling, for complex 
extremals was shown to be a consequence of the Einstein-Hamilton- Jacobi equation for 



complex semiclassical phase 35, 13 . 



4. The homogeneous inflaton mode 



Here we calculate the contribution of the Lorentzian inflaton mode f<^(t) to the 
probability distribution p(Lp, t), given by the preexponential factor Ay^/|t>^(t) | of ( p. 301 ). 
The algorithm (|2.30|) was obtained within the ADM reduced phase-space quantization 
implying that is a physical degree of freedom while the scale factor a and the lapse 



function are determined by the gauge condition ( p.27|) . This gauge condition intro- 
duces cosmic time (with unit lapse A^ = 1) [|1^], therefore the Euclidean and Lorentzian 
wave equations for v^{t) ( ^.31| )-( p.32| ) can be obtained by the same linearization pro- 
cedure as in Sect. 3 (performed in A^ = 1 gauge). After exclusion of 6a of this section. 



this equation boils down to the homogeneous version of Eq. (p.20|) with u^{t) replacing 
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5ip 



w ^ ^ 1-;^-'*^ + o(,4/ki.?)} - (4.1) 



(one should remember that the operator in ( p.20| ) enters the quadratic part of the 
action with the factor AH^ = ^ 27i'^ sin'^ 9 / H) . This equation was, however, 

obtained only in the lowest order approximation while we would need to know f,^(t) in 
a subleading order of the slow roll expansion. The remedy is to use the one-parameter 
family of classical solutions already known in this approximation. The needed linearized 
mode can be obtained by differentiating this solution (|3.24|) with respect to (po 

= (4.2) 

and then analytically continuing it to the Lorentzian spacetime (|2.33|) . Using ( p.21| ) 
one has the Euclidean mode in the first subleading approximation 

ml 1 + 5 f,2 + cos^ 9 , ^ 1 ^^l , 

uJt) = IH i-TT <9 -tan - + 21ncos tan^ - } . (4.3) 

' 2-K^ll-Q^\ l + cos^ 2 2 2 2] ^ ' 

Since it has the form 1 + Oiml/l^lip^ for |^| ^ 1 the modulus of the corresponding 
Lorentzian mode gives the contribution to the exponential of p{ip, t) 

^^T^^=0{mlm^l) (4.4) 

which is by one power of 1/^ smaller in magnitude than the subleading term /i/v^o of 
the corresponding tree- level contribution (|3.43|) . 

Now we have to calculate the Wronskian norm of this mode A^p. From ( [4.1|) it follows 
that the Wronskian operator in the Euclidean time r = 9 / H {W^ =W^— = i W) 
equals 

£1 a9 

The time-independent inner product ( p.34| ) can be calculated at the nucleation point 
t = where it equals 2u^ W^u^. Applying the Wronskian operator to ( [4.3| ) we get 



A 



so that the final contribution of the inflaton mode equals 
Ai/2 

In = const - 21n(^o + 0{ml/\^\^l). (4.7) 
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5. One-loop effective action on the DeSitter instan- 
ton 



In this section we begin calculating the contribution of the one-loop effective action 
( |2.35| ) to the distribution function (p.30| ). As it was discussed in the end of Sect. 2, this 
quantity should be calculated at the extremal Q = Q^°^ + 5Q which differs from the 
exact DeSitter background (with constant inflaton) Q^^^ by the corrections of the slow 
roll expansion: 



ri-ioop= ^Trln '^ ^f^- 



2 6^6^ 



r'l°^ioop + '^r'l.ioop. (5.1) 



So here we obtain F^'^jQ^jp and in the Sect. 6 develop the perturbation theory for (5ri_ioop. 

It is important that in contrast to the 2-dimensional minisuperspace sector (a, (p) 
that was only probed by the tree-level approximation of the theory, the one-loop order 
involves the contribution of all fields inhabiting the model. Without loss of generality 
we shall assume that its low-energy (sub-Planckian) sector is given by the inflaton- 
graviton action ( |2.21|) plus arbitrary set of Higgs scalars x{^)y vector gauge bosons 
Af^{x) and fermions iplx). It can also include gravitino, but we shall mainly focus at 
this sector of spin 0, 1/2, 1 and spin 2 fields. In the full Lagrangian 



-WfIM)-T.^^A +^int(^,X,^M,^) (5.2) 

we single out the interaction of Higgs, vector and spinor fields with the inflaton field 
-^int(v^) X) V^)- ■'■^^ nonderivative part has the form 

^int = ^xV^ + ^S-A^Jv^^ + U^i'tp + derivative coupling (5.3) 
with Higgs A^, vector gauge gA and Yukawa coupling constants. In the Lagrangian 




(5.2) the inflaton field can be regarded as one of the components of one of the Higgs 
multiplets which has a nonvanishing expectation value in the cosmological quantum 
state. In its turn the choice of the interaction Lagrangian here is dictated by the 
renormalizability of the matter field sector of the theory ( |5.2| ) and by the requirement 
of local gauge invariance with respect to arbitrary Yang-Mills group of vector fields A^. 
The terms of derivative coupling in (|5.3| ) should be chosen to guarantee the latter, but 
their form will not be important for the conclusions of this paper. On the contrary. 
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the quantum gravitational effects will crucially depend on the nonderivative part of the 
interaction Lagrangian. 

A very important property of the functional (|0| ) is that it is calculated on shell, 
that is on the solution of classical equations, and therefore, as is well known from the 



theory of gauge fields |^ |38|, |T^, is independent of the choice of gauge conditions 
used for its construction (or ADM reduction to physical variables in terms of which 
it reduces to the functional determinant in the physical sector). This freedom can be 
used to transform ( ^.11) identically to the background covariant gauges in which the 
one-loop action takes the form of the functional determinant of the covariant operator 

F = ^-^14^ (5.4) 
og og 

acting in the full space of gauge and ghost fields 

9 = (fix), x{x), -ipix), Af,{x), V^^(x), g^„y{x), ...,C(x), C{x)). (5.5) 

Here /tot [9] is a total action defined on this space of fields including the covariant 
gauge-breaking and ghost terms, C's denote all possible gauge and coordinate ghosts. 
Spatial components of fields in the nonghost sector of (|5.5|) form the canonical super- 



space of the theory (|2.27|) . As compared to (|2.27|) we only added xi^) ^ the set of all 
scalar multiplets of the model other than the inflaton field (p{x), introduced above. 

The possibility to convert the one-loop action to covariant form is very important, 
because in this form it admits covariant regularization and renormalization of inevitable 
ultraviolet divergences and allows one to obtain correct scaling behaviour, quantum 
anomalies, etc. The full set of gauges for internal gauge symmetries and gravitational 
diffeomorphisms can be chosen in such a way that the operator ( |5.4| ) becomes minimal 
- diagonal in second-order derivatives forming a covariant D'Alambertian □. Moreover, 
on the exact DeSitter background with constant inflaton fleld this operator can as a 
whole be reduced to the block-diagonal form 

F = diag(-n, + X,) (5.6) 

with blocks — + Xg belonging to 0(5) irreducible representations of spin s = 
0, 1/2, 1, 3/2, 2, ... on the 4-dimensional sphere of the Euchdean DeSitter space. For 
constant potential terms Xg of these operators (which is the case of Q*-"-*) their spectra 
are well known [p9|, HOl pi] and can be used for calculating their functional determinants 



under a proper covariant regularization. We shall use a ^-functional regularization p2 
in which a contribution of every block of (|5.6|) equals 



- ^Tr ln(-n. + Xg) = 1C(0) + ^C(O) In (5.7) 
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with the generahzed ^-functions built of dimensionless eigenvalues of the rescaled op- 
erator 



Cs{p) = E ^"^ H-\-n, + X,)Mx) = A0a(x) (5.8) 

(the scale of the DeSitter instanton is its inverse radius H, so that the differential 
operator here is effectively defined on a sphere of unit radius with a correspondingly 
rescaled potential term). As a result the full one-loop action takes the form 



'^'^ -ic'(0) + ^C(0)ln|?. 

C'(o) = E^^C(o), c(o) = E^^C(o). (5.10) 



-raoop = 7TC'(0) + ;,C(0)ln^, (5.9) 



Here the weights Ws - positive and negative integers - reflect the statistics of the field 



and also the details of transition from the original local fields ( |5.5|) to the decomposition 
in irreducible 0(5) representations of (|5.6|) 0. 

In the equations above /i^ is a mass parameter reflecting the renormalization am- 
biguity resulting from the subtraction of logarithmic divergences proportional to C(0)- 
This quantity plays a very important role because it determines the leading high-energy 
behaviour of the one-loop action and correspondingly the anomalous scaling behaviour 
of the distribution function p ~ H~''^^\ As it was observed in [|1^, [1^, ^ it can 
produce a principal quantum cosmological mechanism - to make the distribution func- 
tion of quasi-DeSitter models normalizable in over-Planckian domain and generate the 
inflationary probability peak. Therefore we begin with the calculation of this quantity. 



5.1. The anomalous scaling ({0) 



Apart from the method of 0(5) irreducible representations and Eqs. (^) and 
(|5.10|) the total anomalous scaling ^(0) can be obtained by a more universal Schwinger- 
DeWitt technique 



C(o) 



7], as an integral over the spacetime 



a2[x) 



(5.11) 



of the second Schwinger-DeWitt coefficient a2{x) in the proper time heat kernel expan- 
sion for the operator (|5.4| ). This technique easily allows one to derive the mechanism 



^This transition involves certain Jacobians which have the form of positive and negative powers of 
operators in Q [|9[|^, |l| leading to sign factors of the above type. The example of this procedure 
is presented below for gauge vector field. 
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of suppression of the over-Planckian energy scales due to a big value of the nonmin- 
imal coupling constant |,^| ^ 1. Indeed, for large masses of particles the dominant 
contribution to a2{x) 



«2(x) = - E< + 4E^A-4E^t +••• (5-12) 




is quartic in their masses with the sign factor depending on statistics (and weight 
factors given by the number of the corresponding tensor field components). In the 
model ( |5.2|) -( ^^) on the background with a big constant field (po Higgs scalars, vector 
gauge bosons and fermions acquire by the analogue of the Higgs mechanism the effective 
masses induced by the interaction Lagrangian 



m 



2 



- = g^^Q, = f^if^. (5.13) 



Beeing integrated over the 4-volume of the instanton Stt'^/SH^ these masses generate 
in view of the expression ( ^.5[ ) for H{ipQ) the following dominant contribution to C(0) 

C(0) = Z[1+ 0{ml/\^\^l) ], Z = 6^jA, (5.14) 
where A is a following fundamental combination of matter fields coupling constants [|l| 

^ = ^ f E 4 + 16 E - 16 E 4 V (5-15) 

\ X A i' / 

Thus for large |^| ^ 1 and positive A ~ 0(1) we have Z ^ 1 which is a corner stone 
of the quantum gravitational mechanism that suppresses the over-Planckian energy 
scales, p ~ H^^ ^ 0, H ^ oo and, thus, serves as a justification of the semiclassical 
expansion. 

It is important that this mechanism is entirely generated in the renormalizable 
matter-field sector of the theory consisting of the multiplets of the standard model, 
because the graviton-inflaton sector (|2.21|) yields the contribution independent of ^ 
^ as well as the spin-3/2 gravitino field (considered to be noninteracting 
with inflaton) 

171 589 

C(0)graviton+inflaton 10~' ^^■^S'^^"*™'-' \S0 ' (5.16) 

This property suggests that this mechanism is robust against high energy modifications 
of the fundamental theory designed to solve the problems of nonrenormalizability in 
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perturbative quantum gravity 0. 

Our purpose now is to go beyond the lowest order approximation (|5.14|) and find 
subleading corrections. For this we would need the subleading term in the expression 
for the Hubble constant ( p.5| ) and also use exact expressions for ^(0) of massive scalar, 
vector and spinor fields on the Euclidean DeSitter background. For a scalar field with 
the mass and the constant of nonminimal interaction (which is different from 
this expression reads ^ 

(5.17) 



, , - 29 , 1 1 

Cv(0) = 55-4{, + 12«?--^ + -^. 



Using ( p3| ) and ( [5.13| ) in this equation we have 



Cx(0) 



, 3A2(l + 25)|e| 



m. 



2A. 



A2 ^ 47rA2 
The structure of this expression 



\^\ + 0{ml/vl 



(5.18) 



(5.19) 



demonstrates the nature of the perturbation theory that we shall use in what follows. 
It has two smallness parameters: ^^p/I^Iv^q 1 ^ the parameter of the slow roll 
expansion and ^ 1 - the parameter of this particular model with nonminimal 
inflaton. Below we shall see that |^| ~ 2 x 10^ and the probability maximum in the 
distribution in question will be for mp/ipQ ~ 0.03. This means that for the most 
important range of values of the infiaton field 



1 



(5.20) 



and in the subleading approximation the last two terms in ( |5.18| ) and ( |5.19| ) can be 
discarded. In what follows we shall invariably follow this rule. Thus we have 

1 



Cx(o) = se 



\2 
2 ^ 

A2 



1 



26 ml 



+ 0(1/1^1) 



(5.21) 



4vr lelv^g 

Similar calculations for vector and Dirac spinor fields with masses from ( ^.13 ) give 
the result (taking into account their statistics) 

1 



a(o) 

C^(0) 



2 9 A 



1 + 



2 r 



1 + 



25 ml 

1 + 2(5 ml 
4vr 



+ 0{im) 



(5.22) 



(5.23) 



^Another advantage of this mechanism related to big |^| is that it aUows one to disregard the 
known and thus far unresolved problem of discrepancies between the renormalization in covariant and 
unitary gauges observed on topologically nontrivial curved spaces for fields of spins s > 1 ^ . 
These discrepancies are independent of the nonminimal coupling constant and, thus, negligible for 
lei » 1. 
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Obviously, the contribution of Majorana or Weyl spinor fields is one half of the expres- 
sion (|5:23|) . 

Thus, the total anomalous scaling of the theory on the exact DeSitter background 
reads 



C(0) = 6^A 



l + ii^#. + 0(l/kl) 

471- K|<^0 



(5.24) 



It should be emphasized again that the graviton-inflaton sector does not contribute 
to this expression. We have seen that dependence on ^ in Cx(0) above arises due to the 
terms m^/H^ and m^/H"^. Indeed, in these terms H depends on C,- However, while 
the effective masses m^, m^, of non-inflaton scalar, vector and spinor fields do not 
depend on ^, the effective masses of the graviton and inflaton are ^-dependent (see 
eg. (^.81) for the effective inflaton mass) and this dependence cancels the dependence 
of H on ^. One can show also that the subleading terms in C(0)gravity+infiaton will be 
at most O(ICI) in the terminology of ( ^.19| ) and thus will be discarded in what follows 
From calculations below it will be clear that the same is true for C'(0) foo- Thus in 
the subleading approximation the one-loop part of the probability distribution will be 
contributed by the non-infiaton scalar, vector and spinor fields. 



5.2. Calculating C'(0) 

The calculation of C'(0) requires the knowledge of the finite part of the effective 
action unrelated to its ultraviolet divergences associated with C(0). We shall calculate 
it by the technique of generalized (^-functions developed in EDI, HTI] for fields in the basis 



of irreducible 0(5) representations of different spins s { f).Q\) . In [41| these calculations 
for (^-functions (|5.8| ) were put in a unified framework in the form of a special function 
defined for 3?ep > 2 and analytically continued to p = 

Up) = ^(2s + l)F {p, 2s + l,{s + 1/2)2, ^ (5 25) 

F{pXaM= E ( 2_u,l ^ ^ep>2. (5.26) 
Here the parameters bg are related to potential terms of operators in (|5.6|) according 



to 41 



9 Xo 13 Xi Xi/2 , , 

A~lP'^' = ^~lP' = (5.27) 



Particular values of the function ( |5.26| ) equal 



F(0, k, a, b) = ]b{b - 2a) + ^a{3e + Qk + 2) - ^k\k + 2)^ + (5.28) 
4 24 64 12U 
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k,a,b) = ^b-^- ^k{k + 2)-^{b-a)^(^^ + l±Vb^ 
±y) = ip{x + y) + ip{x — y), 



(5.29) 
(5.30) 



where ip{x) is a logarithmic derivative of Euler's F function. Using the corollary of 

(PD 



^F'{0,k,a,b) = F{l,k,a,b), F' = ^, 
ab dp 

we can find F'{0) and hence C'(0) by integrating ( [5.31| ) 
F'{0,k,a,b) = ^b'-^b-^bk{k + 2) 



1 fk \ 

-- dz{z - a)^ f - + 1 ± v^j + C, 



where the constant C is given by derivatives of Hurwitz functions 



(5.31) 



(5.32) 



C = 2C'h{-^, 1 + k/2) - 2aC^(-l, 1 + fc/2), Ch(x, y) = Y. > ^ w - 

n=o \^ + y) 

The expression ( |5.32|) is rather complicated and generally can be obtained only 
numerically. However, we need only its dependence on (^q contained in b^ and in the 
approximation retaining only the first two terms of the generic expansion ( p.l9| ). Since 
the potential terms of ( p. 6] ) are basically given by effective masses induced according to 



0(v9q), the corresponding parameters b^ ( |5.27| ) have the structure 



&f + ^ + o(i/|ei), bf 



om). 



Using (p:28|)-(|5:29|) and expanding {^1^) and (|53^) in ^i^V^o 

one obtains 

C(0) = C(0)|,(o) + -^-^ + 0{ml/^l), 



C(o) = C(o)L(o) + 



6 (^5 
2s + 1 f 1 r 



(5.33) 



(5.34) 



bf -{s + 1/2)' 



^ s + - ± i 
2 



do) 



m 



+0^-|j. (5.35) 

For the calculation of the full effective action (|5^ ) we would also need ln(/i^/if^). It 
is worth transforming this quantity separately for every irreducible s-component 

,2 

(5.36) 



In^ = ln4 + In^ + ln(-6f) + S- + O 

11 LOr, JVf, 



/i 
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Then, disregarding the (^o-independent part and using in 

the asymptotic expansion for ip{z) at large z |J8|, ip{z) ~ In^; — 1/22; + . . ., we see 



that the logarithmic in b^^^ term of (|5.36|) gets cancelled by the logarithmic term of 



^ (s + 3/2±i( 



1/2 



and the final answer for a partial one-loop action reads 



C(0) + C(0)ln|^ 



const + C(0) In 



+ 



2g + UW6W 
5 V'g 



'^ + lln^ 
4 2 



(5.37) 



Now we can go over to the calculation of the total one-loop action for the model 
.21). For a Higgs scalar field x the contribution reduces to the above equations ( ^.34| ) 

(see eq. ( p.l3| ) for the effective mass of the Higgs 



and (|5.37|) with s = and Xq = 
field), so that in view of the expansion for H ( p.5| 



'Iscal 



3\^{l + 2S) 



m. 



(5.38) 



For a gauge vector field the situation is more complicated because its one-loop action 
in the most convenient gauge V^A^ = (leading to minimal operator with diagonal 
derivatives) equals |^| 



-Tr ln(-n5^; 



m. 



Tr ln(-nc 



- Tr In 
2 



+- Tr ln(- 



Tr ln(-nr 



3H' 



-'0), 



(5.39) 



where the subtracted term is a contribution of ghosts and the other terms represent the 
decomposition of the vector functional determinants into 0(5) irreducible components 
(see footnote after eq. ( |5.1CI| )). Similar procedure holds for a spinor field operator which 
technically must be squared to reduce calculations to that of the ^1/2 + ••• [0- Finally 
we have for vector and spinor fields the coefficients 



do) 



f2 

-12|e|^. 



-12|e| 



3/^(1 + 25) 



3^? 



27rA 
(1 + 26) 



-m 



Pi 



27rA 



m. 



(5.40) 
(5.41) 



and the corresponding Xa, X^ different from the squared masses (|5.13|) by terms that 
go beyond our approximation. Their use in ( ^.34|) and (|5.37| ) allows one to repro- 
duce the expression ( p.24| ) for the total ^(0) obtained above by the Schwinger-DeWitt 
method and get the final algorithm for the total one-loop effective action on the De- 
Sitter background 



r(o) 

l-loop 



const — 3 — A 
A 



1 + 



1 + 25 ml 
4vr 



In 
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where the coefficient A is defined by Eq. (|5.15| ) and S is a following new combination 
of coupling constants 



B = ( E 4 In ^ + 16 E In - 16 E 4 In 4 1 • (5.43) 



6. Perturbation theory for the one- loop effective 
action 

Due to the presence of slow roll corrections 5Q the effective action ( [5.1|) acquires 
the contribution which in the first order of the perturbation theory equals 



(6.1) 



FG{x,x') = 5{x,x') (6.2) 

where G = G{x, x') is the Green's function of the operator F (eq. ( |5.4| )) on a four-sphere 
and 6F is the variation of this operator induced by 6Q = {6a, 6ip). 

The expression ( |6.1| ) is incomplete unless one fixes uniquely the Green's function 
and specifies the functional composition law 6FG in the functional trace. One should 
remember that the kernel G{x,x') is not a smooth function of its arguments and 
its irregularity enhances when it is acted upon by two derivatives contained in 6F. 
Therefore one has to prescribe the way these derivatives act on both arguments of 
G{x, x') and how the coincidence limit of the resulting singular kernel is taken in the 
functional trace 



Tt[6FG] = J d^xtrSF{V,V)G{x,x') 



(6.3) 



(tr denotes the matrix trace operation over tensor indices). The specification of trace 
operation follows from the procedure of calculating the Gaussian path integral over 
quantum disturbances which gives rise to the one-loop functional determinants. As 
was shown in [^] (see also |^0[ and |3^) the functional determinant of the differen- 



tial operator generated by the Gaussian path integral is determined by the variational 
equation 6lnDetF = 6Tt \nF = Tr 6 FG, where the Green's function G{x,x') satis- 
fies the same boundary conditions as the integration variables in the Gaussian integral 
and the functional composition law 6F G implies a symmetric action of spacetime 
derivatives on both arguments of G{x,x'). In what follows we shall describe in detail 
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the calculations for the case of the scalar field and then give the result for other higher 
spin contributions. 

For the field x{^) with the wave operator 

F = -g''''V,V, + ml{^o) (6.4) 



the symmetrized variation of this operator ( |6.4] ) looks like 

SF{V, V) = V^'6g^''{x)V. + 5mJ(<^), (6.5) 
where V^' and Vjy are acting on different arguments of G{x, x') in ( |6.3|) 0. 



The Green's function for a scalar field on a four-sphere is well known , however 



we have to regulate a divergent expression (6^) and will do it by replacing the Green's 



function with its operator power G^^^{x,x') = F ^ ^S{x,x'), p — 0, 

1 

G^+P(x,x') = — -/ dssPexp{-sF)6{x,x'). (6.6) 

T[p + 1) Jo 

The corresponding heat kernel can be constructed by noting that in view of DeSitter 
invariance both the Green's function and its heat kernel are functions of the world 
function a{x, x') - one half of the square of geodetic distance between the points x 
and x' which can also be expressed in terms of the angle 9 = 6{x,x') between these 
points on a sphere of radius R = l/H, a{x,x') = Similarly, the dedensitized 

delta-function above can be constructed in terms of this angular variable y = cos 6 

The scalar D'Alambertian acting on this function of 9 yields the operator of the Leg- 
endre equation after the commutation with the derivative d/dy above 



u — - — 
dy dy 



l-y^)^,-2y— + 2 



d^ ^ d 



(6.8) 



dy"^ dy 

which suggests to expand the delta function in the series of Legendre polynomials Pn{y) 
- eigenfunctions of the Legendre operator with eigenvalues —n{n + 1) 

CO , 1 . 

%-i) = 5:u+- p„(y), (6.9) 

n=0 ^ 



^Note that F enters the action with the metric dependent factor g^^^F, however g^/^ is not 
varied here, because the contribution of this overaU factor is cancehed by the contribution of the 
local measure For this reason, in particular, the effective action is given by the functional 

determinants of F instead of those of g^^'^F. The difference between the corresponding results as well 
as the contribution of local measure are formally proportional to unregulated 5(0) which vanishes in 
dimensional rcgularization, but given by ^(0) in the zeta-functional one. Therefore in the regularization 
we use here these terms require a careful bookkeeping. 
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whence it follows that 

exp{-sF)6{x,x') = f (n + e-^l'^-^'/'^'-'^V^' dP^{y)/dy, (6.10) 

where 60 = 9/4 — rn^R^ coincides with the expression given by ( p.27| ) for a spin-0 case. 
Substituting it into ( |6.6|) , integrating over s and expressing the Legendre polynomials 
in terms of the hypergeometric function we finally obtain 

ri+vi '^ = - (-l)"+Mn + l)(2n + l) 

^ 167r2i?2-2p A. ^ ^/2)2 - 60]'+^ 

1 



xF (^1 - n, n + 2; 2; - + - cos ^2a{x, x')/R^j . (6.11) 

This expression allows one to obtain the coincidence limits of the Green's func- 
tion and its derivatives arising in ( |6.3|) . In view of the well known coincidence limits 
a{x,x) = 0, V/,a(x,x')U'=x = 0, V^V^'a(x, x')U'=x = ~g^lu' we have 

G^-nx X) = y + + (6 12) 

VV,G'^^(rr') 9,u9'/\x) ^ n{n + l){2n + - n){n + 2) 



Infinite series here are calculable by the technique of ^ They have the pole 
structure in p — 0, A/p + B + 0(p), and lead to 5ri_ioop in the form of the finite part 
of the variation ( |6.1| ) []: 

5ri_ioop = \{hA, + /2A2) In 1^ + i(/ifii + /252), (6.14) 
where Ii and /2 are given by the integrals 

and Ai, Bi, i = 1,2, are the following pole and finite parts of the above two sums 



A, = - 6„ (to - i) 



''Note that functional regularization is formally free from divergences - pole terms in p ]4^ . 
These terms are artifacts of the variational equation ( |6.l[ ) which differs from the variation of finite 
C'(0) exactly by the pole term in p. Indeed SCip) = -{1 + pd/dp)TTSF G^+p and for TrJFG^+P = 
A/p + 0{l), p~^0, equals the finite part of Q (5C'(0) = [-Tr SFG^+p + a/p]^^^. 
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Since 6m'^{(p) = X^ipoS(p, Sg^'^ = diag{0,—2g Sa/a), with S(p and 5a given by 
eqs. (|3.21|) and ( |3.22|) , the final resuh for a scalar field contribution to 5ri„ioop reads 
up to terms 0{mp/(fl) 



loop 



<«)) 



5*C.(0)ln:^, 



ttAVo 



A 



vr TT In 2 1 

K = — + ^ 0.2. 

96 12 72 



(6.18) 
(6.19) 



The leading terms of S({0) and S('{0) for vector and spinor fields have a similar 
form and when composed with the contribution of Higgs multiplets above give rise to 
the total (5ri_ioop- Similarly to the unperturbed part (|5.42|) it expresses as a function 
of the universal combinations of coupling constants A and B given by Eqs.( ^.15| ) and 

(Pi 



5Tu 



loop 



kB. 



(6.20) 



7. Probability maximum of the distribution func- 
tion 



Combining the tree-level part ( |3.43| )-( p.45[ ) with the contributions of the infiaton 
mode ( [4.7| ) and perturbative contributions of the one-loop effective action ( |5.42| ) and 
(|6.20|) we finally arrive at the distribution function ( |2.30 ) of infiationary cosmologies 



PT,NB{fo) = N exp 



m 

A 



m. 



m. 



^0 W^m, 



(7.1) 



where is a field independent normalization factor and a± and (3 are the following 
functions of coupling constants of the model 

6 



±^aA + B)±3.B'— 



1 + 26 ,1 + 5 

— 6K . 



An 71 

involving the parameter 6 ( p.23|) and universal combinations A and B. 



(7.2) 
(7.3) 
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The second term in the exponential of ( [7. 1|) confirms the conclusions of [|ri|, [12], |l| 
that loop corrections can drastically change the predictions of classical theory and sup- 
press the contribution of the over-Planckian energy scales due to big positive anomalous 
scaling of the theory 

Pt,nb{vo) ~ V^o V^o oo, (7.4) 



with Z given by eq. ( [5.14 ). Extra power —2 of (fo here comes from the contribution of the 



inflaton mode ([4.7|) neglected in within the 0{w?p/ip1) accuracy but important for 
LfQ oo. For positive Z = 6|^pA/A this asymptotics can be regarded as a justification 
of a semiclassical expansion. 

The equation for the extremum of the obtained distribution at (po = (fj can be 
represented in the form 




where plus or minus signs correspond to the tunnelling or no-boundary wavefunctions 
respectively. To analyze the existence of its solution we shall have to prescribe a certain 
reasonable range of parameters and try solving it by iterations. In the next section 
we shall briefly discuss the model of nonminimal inflation with big |^| ~ 2 x 10^ ^ 
1 and S = 0(1) that was used in our previous work [0 as a good candidate for a 
quantum origin of the early Universe at a sub-Planckian energy scale (around GUT 
scale). In we only took into account the tree-level part and the leading logarithmic 
behaviour of the one-loop effective action, which correspond to retaining in ( [7.1D only 
the first two terms of the exponential with the parameter a truncated to the first {A 
and independent) term of ( [7.2|) . In it was shown that the requirement of the 



minimal admissible duration of the inflationary stage (|2.5|) imposes upper bound on 
the combination of coupling constants A 1.3. Here we shall show that the qualitative 
estimates of ||l|] remain also true after the inclusion of perturbative corrections obtained 
above. 

The justification of the results of consists in the observation that the second 
term in the right-hand side of ( |7.5| ) can be regarded small and taken by perturbations. 
Indeed, from the upper bound on A we can assume that both A and B are of the 
order of magnitude one 

A = 0(1), B = 0{1). (7.6) 



This follows from the comparison of expressions ( ^.15|) and (^.43|) and a natural as- 



sumption that the estimate for B follows from that for A, unless strong cancellation 
takes place between different separately big terms of the expression ( |5.15| ) for A. From 
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these bounds we see that a± = 87r(l + 6) + 0(1/10), /? = 0(1/10) due to the numerical 
values of the coefficients in (|7.2| )- (|7.3| ) and the fact that 6 = 0(1). Then it follows that 
l3A/a± ~ A/80'7r(l + 5) ~ A/250 and, therefore, the second term in the right-hand 
side of ( [775] ) can indeed be treated perturbatively with a good smallness parameter 
A/250. With the same precision the expression for a± reduces in the leading order to 
a = 87r(l + (5). It is needless to say that for smaller values of A and B this approxima- 
tion works even better, because the leading term in the right-hand side of (|7.5| ) grows 
and the corrections decrease for A, B ^ 0. 

For 6 > —1 (which is the case of a standard classical scenario with a finite duration 
of inflation, see Sect. 2) the leading order solution of ( [7.5| ) exists only for the case of 
a tunneling wavefunction (sign plus) and coincides with the result of |ll|. Thus, the 
parameters of the inflation probability maximum - the mean value of the inflaton field 

and its quantum dispersion A^p = [—d'^ln p{Lpi) / dLpf\~^^'^ are 



For the no-boundary quantum state of the Universe the peak can be realized only for 
6 < —1 and, thus corresponds to the classical scenario with endless inflation stage ]l|]. 



8. Nonminimal inflation and particle physics of the 
early Universe 

The present state of inflation theory is consistent with observations of the cosmic 
microwave background radiation anisotropy in the COBE [|14| and Relikt satellite 
experiments. In the chaotic inflationary model with a nonminimal inflaton field ( p.21| ) 
the spectrum of perturbations compatible with these measurements can be acquired 
in the range of coupling constants A/^^ ~ 10~^° [|1^, |5^ (the experimental bound 
on the gauge-invariant density perturbation P,^(A;) = A''|(A/^^)/87r^ in the fc-th 
mode "crossing" the horizon at the moment of the e-foldings number N^). The main 
advantage of this model is that it allows one to avoid an unnaturally small value of A in 
the minimal inflaton model [@] and replace it with the GUT compatible value A ~ 0.05, 
provided ^ ~ —2x10^ is chosen to be related to the ratio of the Planck scale to a typical 
GUT scale, |,^| ~ nip/v. For these values of coupling constants the parameter ( p.23| ) 
is 5 ~ Snv/nip ~ 10"'^ (thus easily satisfying a much weaker upper bound 6 = 0(1) 
assumed above). As far as it concerns the Hubble parameter (|3.5| ) and the number of 
e-foldings ( p.24| ) at the inflationary peak ( |7.7| ) obtained above, in the leading order in 
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1^1 S> 1 they are given by 



^(¥^/)=mp^^— (8.1) 

N{^,) = ^ (8.2) 

and satisfy the bound N{(fi) > 60 with a single restriction on A, A < 1.3. This 
restriction gives rise to the bounds ( |7.6|) and, therefore, justifies the smallness of the 
perturbative corrections of the above type, that were neglected in our previous paper 
[|I|. These perturbative corrections contained in the second term of the equation ([7^), 
unfortunately, violate the conclusion of [jT| on the absence of renormalization ambiguity 
in the energy scale of inflation, but as we see this ambiguity in the choice of /^^ is 
strongly suppressed by the smallness of the ratio (3 A/ a. 

On the other hand, the restriction on A justifles a slow-roll approximation, because 
the corresponding smallness parameter is ^/Hip —A/dGn'^ ~ — 10~^. For the above 
small value of 5 <^ 1 and A ~ 1, the obtained numerical parameters describe extremely 
sharp inflationary peak at ipi with small width and sub-Planckian Hubble constant 

(fj ~ 0.03mp, a ~ lO^W, H{ifj) ~ 10"W, (8.3) 

which is the most realistic range of the inflationary scenario. The smallness of the width 
does not, however, lead to its quick quantum spreading: the commutator relations for 
operators and (p, [0,0] — 2/(127r^|^|a^) [0, give rise at the beginning of inflation, 
a ~ H^^, to a neglig ible dispersion of 0, A0 ~ /127r'^\^\a ^ {8/ A){y/\/\^\)\0\ < 
\0\. It is remarkable that the relative width 

Aif AH . , , 

-T. ^ __ ^ 10-5 8.4 

Lfj H 

corresponds to the observable level of density perturbations, although it is not clear 
whether this quantum dispersion Ayj is directly measurable now, because of the stochas- 
tic noise of the same order of magnitude generated during the inflation and superim- 
posed upon Aip. 

All these conclusions are universal for a generic low-energy model ( [5.2| ) and (apart 
from the choice of |^| and A) universally depend on one parameter A of the particle 
physics model. This quantity should satisfy the lower and upper bounds 

^ < A < 1.3 (8.5) 

in order respectively to render Z positive, thus suppressing over-Planckian energy 
scales, and provide sufficient amount of inflation (A should not, certainly, be exceed- 
ingly close to zero, not to suppress the dominant contribution of large |^| in ( |5.14| )). 
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In [H the conclusion was made that this bound suggests the quasi-supersymmetric 
nature of the particle model, because supersymmetry can constrain the values of the 
Higgs A^, vector gauge qa and Yukawa couplings so as to provide a subtle balance 
between the contributions of bosons and fermions in ( |5.15| ) and fit the quantity A into 
a narrow range ( |8.5| ). Now, however, with the inclusion of corrections that go beyond 
the estimates of we have also to provide the boundedness of B ( |5.43| ) which is 
less obvious to be compatible with the bound ( ^.5|) , unless all the terms of ( ^.15| ) and 
(|5.43|) are separately small due to small values of all coupling constants. In the latter 
case the supersymmetry is not needed to explain the restrictions (|8.5|) on the choice 
of a particle physics model. Still, supersymmetry remains a reasonable conjecture 
consistent with this selection criterion ( p.5| ) and sounds coherent with conclusions of 
[p5| where supersymmetry was argued in the opposite case of small |^|. 



9. Conclusions 

Thus, the same mechanism that suppresses the over-Planckian energy scales also 
generates a narrow probability peak in the distribution of tunnelling inflationary uni- 
verses and is likely to suggest the (quasi)supersymmetric nature of their particle con- 
tent. It seems to be consistent with microwave background observations within the 
model with a strongly coupled nonminimal infiaton field. A remarkable feature of this 
result is that it is mainly based on one small parameter - the dimensionless ratio of 
two major energy scales, the GUT and Planck ones, given by the combination of the 
coupling constants \/~\/\^\ ~ 10~^. 

Big value of |^| is actually responsible for the fact that the one- loop corrections 
qualitatively change the tree-level behaviour and produce the inflationary peak of the 
above type. In the absence of powerful nonperturbative methods there is no rigor- 
ous proof that the inclusion of multi-loop orders will not destroy this one- loop effect. 
However, the following qualitative arguments support the conjecture that this will not 
happen. Point is that the effective gravitational constant in this model is inverse pro- 
portional to rrip + 87r|^|(y9^ and, thus, large |^| might improve the loop expansion ||2l|] by 
suppressing the contribution of multi-loop orders. On the other hand, the power-law 
mechanism (|7.4|) of suppressing the high-energy scales is independent of the renormal- 
ization ambiguity (the parameter /i^). This also gives a hope that the obtained effect 
is robust against inclusion of multi-loop corrections. All these corrections are weighted 
by this suppressing one-loop factor and are small at least in some range of (p above the 
probability peak value 99/ where all their curvature invariants have GUT values rather 
than Planckian ones. Therefore, irrespective of what happens at Planckian scales this 
peak at GUT scale will be separated from the unknown Planckian domain by the re- 
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gion with very small values of the distribution function and, at least heuristically, will 
not interfere with nonperturbative quantum gravity or its more fundamental (stringy) 
generalization. 

Obviously, the large value of |^| at sub-Planckian (GUT) scale requires explanation 
which might be based on the renormalization group approach (and its extension to 
non-renormalizable theories [^). As shown in [^, quantum gravity with nonminimal 
scalar field has an asymptotically free conformally invariant (^ = 1/6) phase at over- 
Planckian regime, which is unstable at lower energies. It is plausible to conjecture that 
this instability can lead (via composite states of the scalar field) to the inversion of the 
sign of running and its growth at the GUT scale, thus making possible the proposed 
inflation applications. 

From the viewpoint of the theory of the early universe, the obtained results give a 
strong preference to the tunnehng quantum state. Debates on advantages of the tun- 
neling versus no-boundary wavefunction has a long history ^ ^ |^. At present, 
in the cosmological context the tunnelling proposal seems to be more useful and con- 
ceptually clearer than the no-boundary one, because for its interpretation one should 
not incorporate vague ideas of the third quantization of gravity which inevitably arise 
in the no-boundary case: splitting the Lorentzian wavefunction in positive and nega- 
tive frequency parts and separately calculating their probability distributions. On the 
other hand, the formulation of the tunneling proposal is not so aesthetically closed, for 
it involves imposing outgoing condition after the potential barrier, the unit normal- 
ization condition - before the barrier at a = 0, the requirement of the normalizability 
in variables / |Q, etc. And all this in contradistinction with the closed path- integral 
formulation of the no-boundary proposal, automatically providing many of the above 
properties. On the other hand, outside of the cosmological framework, in particular, 
within the scope of the wormhole and black hole physics, the tunneling proposal seems 
to be helpless. Moreover, at the overlap of the cosmological framework with the theory 
of the virtual black holes it leads to contradictions signifying that the quantum birth 
of bigger black holes is more probable than of the small (Planckian) ones [^. All 



these arguments can hardly be conclusive, because it might as well happen that the 
difference between the no-boundary and tunnelling wavefunctions should be ascribed 
to the open problem of the correct quantization of the conformal mode. Note that 
the normalizability criterion for the distribution function and its algorithm ( ^.30] ) do 
not extend to the low-energy limit (f ^ 0, where a naively computed no-boundary 
distribution function blows up to infinity, the slow-roll approximation becomes invalid, 
etc. This is a domain related to a highly speculative (but, probably, inevitable) third 
quantization of gravity [^], which goes beyond the scope of this paper. Fortunately, 
this domain is also separated from the obtained inflationary peak by a vast desert 
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with practically zero density of the quantum distribution, which apparently justifies 
our conclusions disregarding the ultra-infrared physics of the Coleman theory of baby 
universes and cosmological constant . 
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Figure captions 



Fig.l Graphical representation of the Lorentzian spacetime L nucleating at the bounce 
surface from the Euclidean manifold E of the no-boundary type, having the topol- 
ogy of the four-dimensional ball. 

Fig. 2 Two-dimensional minisuperspace of the scale factor a and the inflaton field in 
the chaotic inflation model. The Euclidean extremal (in the slow roll approximation) 
starts at a = with large initial value do in the form of a trajectory that is reflected 
from the caustic S^, a ~ 1/H{4>), and enters the region — > oo, a — > 0. 

Fig. 3 Graphical representation of calculating the quantum distribution of tunnelling 
Lorentzian universes: a composition of the combined Euclidean-Lorentzian spacetime 
iW_ U L with its orientation reversed and complex conjugated copy U L* results 
in the doubled Euclidean manifold 2M - the gravitational instanton carrying the Eu- 
clidean effective action of the theory. The cancellation of the Lorentzian domains L and 
L* reflects unitarity of the theory in the physical spacetime of Lorentzian signature. 
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